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Abstract 

Many extensions of the Standard Model posit the existence of heavy particles with long lifetimes. 
In this Letter, results are presented of a search for such particles, which decay at a significant distance 
from their production point, using a final state containing charged hadrons and an associated muon. 
This analysis uses a data sample of proton-proton collisions at = 7 TeV corresponding to an 
integrated luminosity of 4.4 fb~^ collected in 201 1 by the ATLAS detector operating at the Large 
Hadron Collider. Results are interpreted in the context of i?-parity violating supersymmetric scenarios. 
No events in the signal region are observed and limits are set on the production cross section for 
supersymmetric particles, multiplied by the square of the branching ratio for a neutralino to decay to 
charged hadrons and a muon, as a function of the neutralino lifetime. To allow these limits to be used 
in a variety of models, they are presented for a range of squark and neutralino masses. 
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1. Introduction 

Several extensions to the Standard Model predict the 
production of heavy particles with lifetimes that may be 
of order picoseconds up to about a nanosecond at 
the Large Hadron Collider (LHC). One such scenario is 
gravity-mediated supersymmetry (SUGRA [5]) with R- 
parity violation (RPV) [31 0]. If the results of i?-parity 
conserving supersymmetry searches at the LHC (see, for 
example, Refs. E]) continue to disfavour light super- 
partners, then scenarios in which i?-parity is violated may 
be required if supersymmetry is the solution to the hi- 
erarchy problem [7]. The present (largely indirect) con- 
straints on RPV couplings T would allow the decay of 
the lightest supersymmetric particle as it traverses a par- 
ticle detector at the LHC. Signatures of heavy, long-lived 
particles also feature in models of gauge-mediated super- 
symmetry breaking [B], split-supersymmetry ^9j, hidden- 
valley [TU], dark-sector gauge bosons [TT] and stealth su- 
persymmetry [12j though some of these models are dis- 
favoured [13] by the recent results on the search for the 
Higgs boson [Ml [15] . 

This Letter presents the results of a search for the de- 
cay of a heavy particle, producing a multi-track vertex 
that contains a muon with high transverse momentum (pt) 
at a distance between millimeters and tens of centimeters 
from the pp interaction point. As in the earlier ATLAS 
work |16j . the results are interpreted in the context of an 
i?-parity violating supersymmetric scenario. In this model, 
the signature under consideration corresponds to the de- 



cay of the lightest supersymmetric particle, resulting in a 
muon and many high-px charged tracks originating from 
a single displaced vertex (DV). This can arise from a dia- 
gram such as that shown in Fig.[l] where the decay occurs 
due to the non-zero RPV coupling Ajj^ . 

Searches for related signatures have also been made at 
the Tevatron with y/s = 1.96 TeV pp coUisions. The D0 
collaboration has searched for a long-lived neutral particle 
decaying into a final state containing two muons T7j or 
a bb pair |18j . Similar searches have also been made at 
LEP [19] . The result reported here improves on the previ- 
ous ATLAS search for this signature [HT in several ways. 
The dataset is around 100 times larger, corresponding to 
an integrated luminosity of 4.4 fb~^ [2D]. Furthermore, 
the signal efficiency is increased by improving the recon- 
struction efficiency for reconstructing tracks that do not 
originate from the primary vertex. The dependence on 
simulation to estimate background levels is also removed. 



2. The ATLAS detector 

The ATLAS detector [2l] is a multipurpose apparatus 
at the LHC. The detector consists of several layers of sub- 
detectors. From the interaction point (IP) outwards there 
is an inner detector (ID), measuring the tracks of charged 
particles, electromagnetic and hadronic calorimeters, and 
an outer muon spectrometer (MS). 

The ID is immersed in a 2 T axial magnetic field, and 
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Figure 1: Diagram of the massive lightest neutralino decaying 
into a muon and two jets via a virtual smuon, with lepton-number 
and i?-parity violating coupling X'^^^ . 

extends from a radiu^of about 45 mm to 1100 mm and 
to \z\ of about 3100 mm. It provides tracking and vertex 
information for charged particles within the pseudorapid- 
ity region \ri\ < 2.5. At small radii, silicon pixel layers and 
stereo pairs of silicon microstrip detectors provide high res- 
olution pattern recognition. The pixel system consists of 
three barrel layers, and three forward disks on either side 
of the interaction point (the beam pipe envelope is at a ra- 
dius of about 30 mm). The barrel pixel layers, which are 
positioned at radii of 50.5 mm, 88.5 mm, and 122.5 mm are 
of particular relevance to this work. The silicon microstrip 
tracker (SCT) comprises four double layers in the barrel, 
and nine forward disks on either side. A further tracking 
system, a transition-radiation tracker (TRT), is positioned 
at larger radii. This device is made of straw-tube elements 
interleaved with transition-radiation material; it provides 
coverage within the region |?7| < 2.0. 

The calorimeter provides coverage over the pseudora- 
pidity range |r/| < 4.9. It consists of a lead/liquid-argon 
electromagnetic calorimeter, a hadronic calorimeter com- 
prising a steel and scintillator-tile system in the barrel re- 
gion and a liquid-argon system with copper and tungsten 
absorbers in the end-caps. 

Muon identification and momentum measurement is 
provided by the MS. This device has a coverage in pseu- 
dorapidity of |?7| < 2.7 and is a three-layer system of gas- 
filled precision-tracking chambers. The pseudorapidity re- 
gion \ri\ < 2.4 is additionally covered by separate trigger 
chambers, used by the hardware trigger for the first level 
of triggering (level-1). The MS is immersed in a magnetic 
field which is produced by a set of toroid magnets, one 
for the barrel and one each for the two end-caps. Online 
event selection is performed with a three-level trigger sys- 
tem. It comprises a hardware-based level-1 trigger, which 
uses information from the MS trigger chambers and the 



^ ATLAS uses a right-handed coordinate system with its origin at 
the nominal IP in the centre of the detector and the z-axis along the 
beam pipe. The x-axis points from the IP to the centre of the LHC 
ring, and the j/-axis points upward. Cylindrical coordinates (r, (f>) 
are used in the transverse plane, being the azimuthal angle around 
the beam pipe. The pseudorapidity is defined in terms of the polar 
angle 9 as rj = — lntan(6/2). 



calorimeters, followed by two software-based trigger lev- 
els. 

3. Data and simulation 

The data used in this analysis were collected between 
March and October, 2011. After the application of beam, 
detector, and data-quality requirements, the integrated lu- 
minosity considered corresponds to 4.4 fb^^, with an un- 
certainty of ±3.9% gU]. 

Signal Monte Carlo (MC) events are generated with 
Pythia 6 [12], using the MRST LO** [23, set of parton 
distribution functions (PDFs). Processes are simulated in 
which a qq, q q, or qq pair is produced in the pp collision, 
with each squark (antisquark) decaying into a long-lived 
lightest neutralino and a quark (antiquark). Degeneracy 
of the first and second generations and for the left-handed 
and right-handed squarks is assumed. The masses of the 
gluino, sleptons and third-generation squarks are set at 
such a high value (5 TeV) that they are not directly pro- 
duced in the supcrsymmetry scenario considered here. 

The parameter settings for the samples of signal MC 
events are summarised in Table [T] The chosen values 
of squark and neutralino masses correspond to a wide 
range in the quantities to which the signal efficiency is 
most sensitive: neutralino speed and final-state multiplic- 
ities (see Section [5]) . The Higgsino-gaugino mixing deter- 
mines the neutralino mass values. The signal MC samples 
are labelled in Table [l] according to the squark mass and 
neutralino mass, respectively: MH (medium-mass squark, 
heavy neutralino), ML (medium-mass squark, light neu- 
tralino), and HH (heavy squark and heavy neutralino). It 
is assumed that all RPV couplings other than are 
zero. This leads to each neutralino decaying to fi~ud (or 
the charge-conjugate state). Decays of neutralinos to a 
neutrino and jets final state, which would also occur via a 
non-zero A211, are not considered in this work. The values 
of the X211 coupling in the various samples are selected 
such that a significant fraction of neutralino decays occur 
in the detector volume considered in this analysis. The 
average neutralino proper flight distances {ctmc, where 
tmc is the generated neutralino lifetime) are provided in 
Table HI 

The PROSPINO [H] program is used to calculate the 
signal cross section (hereafter referred to as the supersym- 
metry production cross section) at ncxt-to-leading order; 
cross section values are listed in Table [l] The CTEQ6.6 
PDF set 125] is used. Owing to differences in the pre- 
dicted cross sections for the qq and qq processes from 
Pythia and PROSPINO, and the possible sensitivity 
of the neutralino-speed distribution to the assumed pro- 
cess, signal samples are generated for each process sepa- 
rately with Pythia and reweighted according to the rela- 
tive cross sections of the two processes as estimated with 
PROSPINO. The reweighting procedure, however, has a 
minimal effect on the distributions of physical quantities 



2 



Sample 


rriq 


(T 




(7/3),o 


ctmc 






[GeV] 


[fb] 


[GcV] 




[mm] 


xlO" 


MH 


700 


66.4 


494 


1.0 


78 


0.3 


ML 


700 


66.4 


108 


3.1 


101 


1.5 


HH 


1500 


0.2 


494 


1.9 


82 


1.5 



one SCT hit that is expected but not found i.e. in an ac- 
tive detector element with no known read-out problems. 
Furthermore, the track must satisfy an |?7|-dependent re- 
quirement on the number of TRT hits. No pixel-hit re- 
quirements are applied to the muon track. The combina- 
tion of requirements described above is referred to as the 
muon-selection criteria. 

In order to reconstruct a DV it is first necessary to se- 
lect high quality tracks. This is done by requiring that 
candidate tracks have two or more associated SCT hits 
and a value of \do\ greater than 2 mm. Studies made 
with the MC simulation show that 98% of all tracks orig- 
inating from the primary pp interaction are rejected by 
the selection on |do|. Standard ATLAS tracking [21] is 
based on three passes in which the initial track candidates 
are formed in different ways: initially found in the silicon 
detectors, initially found in the TRT detector, and only 
found in the TRT detector. The algorithms all assume 
that tracks originate from close to the PV, and hence have 
reduced reconstruction efficiency for signal tracks which 
originate at a DV. To counter this problem and recover 
some of these lost tracks, the silicon-seeded tracking al- 
gorithm is re-run with looser requirements on the radial 
and z impact parameters, and on the number of detector 
hits that can be shared among more than one seed track. 
To reduce the rate of false seed tracks, it is required that 
these additional tracks have > 1 GeV, which is greater 
than the standard-tracking requirement of pt > 400 MeV. 
This procedure is termed "re-tracking" . Fig. [2] shows the 
efficiency for vertex reconstruction (described later) as a 
function of the radial position of the vertex r^v when us- 
ing standard tracking and re-tracking for sample MH. As 
can be seen, there is a substantial improvement due to re- 
tracking at values of tdv greater than about 10 mm. The 
dips in the plot correspond to losses in efficiency for de- 
cays immediately before a pixel layer, where many tracks 
from the vertex have shared pixel hits and therefore fail 
the selection. 
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Figure 2: The vertex reconstruction efficiency as a function of tdv 
for the MH sample, with and without re-tracking. All selection cri- 
teria, except for the material veto, have been applied. 



Table 1: The parameter values for the three signal MC samples 
used in this work: the assumed squark mass, production cross section 
(calculated with PROSPINO), neutralino mass, average value of the 
Lorentz boost factor (from Pythia), average proper decay length, 
and value of ^211 coupling used in the generation of the signal MC 
samples. 



used in this work and typically leads to changes of much 
less than 1%. 

Each generated event in the signal samples is pro- 
cessed with the GEANT4-based [26] ATLAS detector sim- 
ulation [27] and treated in the same way as the collision 
data. The samples include a realistic modelling of the ef- 
fects of multiple pp collision per bunch crossing (pile-up) 
observed in the data. 

4. Vertex reconstruction and event selection 

It is required that the trigger identifies a muon candidate 
with transverse momentum > 40 GeV and \r]\ < 1.05, 
i.e., in the barrel region of the MS. The latter selection 
is necessary since there is no ID track requirement in the 
trigger and there is a higher fake-muon rate in the (larger 
I77I) end-cap region. 

Owing to the higher instantaneous luminosity of the 
LHC in 2011 compared to that in 2010, most events 
contain more than one primary vertex (PV). The PV 
with the highest sum of p^ of the tracks associated to 
it is required to have at least five tracks and a z posi- 
tion in the range Izpvl < 200 mm. In order to reduce 
cosmic-ray background, events are rejected if they con- 
tain two muons which appear back-to-back. A selection 
\/{tt - (01 - 02))^ - im + V2Y) > 0.1 is applied where 
4'ii 4*2, Vi ^i^d ?72 are the azimuthal angles and pseudo- 
rapidities of the two reconstructed muons. 

A muon candidate is required to have been recon- 
structed in both the MS and the ID with transverse mo- 
mentum pt > 50 GeV (which is well into the region 
where the trigger efficiency is approximately independent 
of the muon momentum) and \r]\ < 1.07. The impact 
parameter relative to the transverse position of the pri- 
mary vertex (cZq) is required to satisfy |dol > 1-5 mm. 
To ensure that the muon candidate is associated with the 
muon that satisfied the trigger requirement, the selection 
•\/ (A(/))2 -I- (A?7)2 < 0.15 is imposed, where A0 (A77) is the 
difference between the azimuthal angle (pseudorapidity) of 
the reconstructed muon and that of the muon identified by 
the trigger. The ID track associated with the muon candi- 
date is required to have at least six SCT hits, and at most 
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Using an algorithm based on the incompatibility-graph 
approach [55], DVs are sought with the selected tracks. 
The method adopted is similar to that used in Ref. 
The algorithm starts by finding two-track seed vertices 
from all pairs of tracks; those that have a vertex-fit 
of less than 5.0 per degree of freedom are retained. A 
seed vertex is rejected if at least one of its tracks has hits 
between the vertex and the PV. Multi-track vertices are 
formed from combinations of seed vertices in an iterative 
process. The following method is used to do this. If a track 
is assigned to two different vertices, the action taken de- 
pends on the distance D between the vertices, li D < Sao, 
where cr is the estimated uncertainty on D, a single vertex 
is formed from the tracks of the two vertices. Otherwise, 
a track is associated with that vertex for which the track 
has the smaller value. In the event of the of a track 
relative to the resulting vertex exceeding 3.0 per degree 
of freedom, the track is removed from the vertex, and the 
vertex is refitted. This process continues until no track is 
associated with more than one vertex. Finally, vertices are 
combined and refitted if they are separated by less than 
1 mm. 

To ensure a good quality DV fit, the x^ per degree of 
freedom of the fit must be less than 5. Furthermore, the 
DV position is required to be in the fiducial region, de- 
fined as I^Dvl < 300 mm and rov < 180 mm, where 
tdv and zdv a-re the radial and longitudinal vertex po- 
sitions with respect to the origin. To minimise back- 
ground coming from the PVs, the transverse distance 
\/ (xdv — xpv)"^ + (j/DV — VwY between any of the PVs 
and the DV must be at least 4 mm. Here x and y are 
the transverse coordinates of a given vertex, with the sub- 
scripts PV and DV denoting the type of vertex. The num- 
ber of tracks N^y associated with the DV is required to be 
at least five. This reduces background from random com- 
binations of tracks and from material interactions. Back- 
ground due to particle interactions with material is sup- 
pressed by requiring mnv > 10 GeV. Here, mnv is the 
invariant mass of the set of tracks associated with the 
DV, using the charged-pion mass hypothesis for each track. 
Candidate vertices which pass (fail) the mov > 10 GeV 
requirement are hereafter referred to as being high-mDV 
(low-mDv) vertices. 

The typical position resolution of the DV in the signal 
MC samples is tens of microns for r^v and about 200 /xm 
for Zdv near the interaction point. For vertices beyond 
the outermost pixel layer, which is located at a radius of 
122.5 mm, the typical resolution is several hundred mi- 
crons for both coordinates. 

Low-mov vertices from particle-material interactions 
are abundant in regions in which the detector material is 
dense. High-mDv background may arise from the random 
spatial coincidence of a material-interaction vertex with a 
high-pT track, especially when this track and the parti- 
cle that initiated the material-interaction vertex originate 
from different primary interactions, resulting in a large 
angle between their momentum vectors. 



To reduce this source of background, vertices that are 
reconstructed in regions of high-density material are re- 
moved. The high-density material was mapped using low- 
TTiDV material-interaction candidate vertices in data and 
true material-interaction vertices in minimum-bias MC 
events. The zdv and tdv positions of these vertices are 
used to make a two-dimensional material-density distribu- 
tion with a bin size of 4 mm in zpy and 1 mm in tdv- 
It has been demonstrated [H] that the detector simula- 
tion describes the positions of pixel layers and associated 
material reasonably well, while the simulated position of 
the beampipe is shifted with respect to the actual posi- 
tion. The use of data events to construct the material 
map therefore ensures that the beampipe material is cor- 
rectly mapped, while the use of the simulation provides a 
high granularity of the map at the outer pixel layers, where 
material-interaction vertices in the data are comparatively 
rare. Material-map bins with vertex density greater than 
an tdv" and ZDV-dependent density criterion are charac- 
terised as high-density-material regions. These make up 
34% of the volume at |zdv| < 300 mm, tdv < 180 mm. 

To ensure that the muon candidate is associated with 
the reconstructed DV, the distance of closest approach of 
the muon with respect to the DV is required to be less than 
0.5 mm. The aforementioned selections are collectively 
referred to as the vertex-selection criteria. 

5. Signal efficiency 

The signal efficiency depends strongly on the efficien- 
cies for track reconstruction and track selection, which are 
affected by several factors: (1) The number of tracks orig- 
inating from the DV and their total invariant mass in- 
crease with the neutralino mass. (2) More tracks fail the 
I do I > 2 mm requirement for small rov or if a neutralino 
with large has a small decay length, thereby leaving 
its daughters pointing back closer to the PV. (3) The effi- 
ciency for reconstructing tracks decreases with increasing 
values of |c?o|. Because the MH and HH samples have the 
same neutralino mass, but different boosts, a cross-check 
of the efficiency estimation procedure can be made. The 
HH sample is reweighted vertex-by-vertex such that the 
versus r] distributions match the MH sample. It was 
checked that they have the same efficiency as a function 
of decay position. 

The total efficiency for each of the signal MC samples is 
shown in Fig. [3]as a function of cr. This is an overall event 
efficiency, derived from the determination of (the single) 
vertex efficiency, and is based on having two displaced ver- 
tices per event in the signal MC. This efficiency as a func- 
tion of cr is calculated from the samples generated with a 
single-cT value (see Table [T]) using a two-dimensional map 
in (I^Dvli ^Dv) (Fig-|4]) and distributions of decay positions 
for different values of cr. Events are reweighted according 
to their probability of being found at different positions in 
the map, where the probability depends on cr. Figure [3] 
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includes systematic corrections and uncertainties that are 
discussed in Section [3 
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Figure 3: The event selection efficiency as a function of ct for the 
three signal samples. The total uncertainties on the efficiency are 
shown as bands; the estimation of the uncertainties is discussed in 
Section El 



ATLAS simulation \^ = 7 TeV 




Figure 4: The efficiency as a function of r^v ^md z^y for vertices 
in the signal MC sample MH. The blank areas represent regions of 
dense material which are not considered when looking for DVs. Bin- 
to-bin fluctuations due to limited MC sample size are compensated 
for by recalculating this map many times, varying each bin content 
within its statistical uncertainty. 



6. Background estimation 

Spurious high-rriDv? high- multiplicity vertices in non- 
material regions could come from one of two sources: (1) 
Purely random combinations of tracks (real or fake) . This 
type of background is expected to form the largest con- 
tribution at small radii where the track density is high- 
est. (2) Vertices from real hadronic interactions with gas 
molecules. Although most of these events will have masses 



below the 10 GeV requirement, the high-mass tail of the 
distribution indicates a potential background, in particu- 
lar if the vertex is crossed by a random track (real or fake) 
at large crossing angle. 

The first source of background mentioned above gives 
rise to displaced vertices inside the beampipe (where there 
are few real hadronic interactions due to the good vacuum, 
but where random combinations are expected to domi- 
nate). The expected background yield is evaluated using 
two independent control samples. One control sample uses 
jet-triggered events (rather than muon-triggered events), 
and the other uses vertices with 4 GeV < mnv < 10 GeV, 
which is below the standard vertex selection requirement. 
This procedure gives an estimate of 0.001^q'°^ vertices for 
this source of background. 

To estimate the second source of background i.e. ver- 
tices arising from real hadronic interactions, including 
those crossed by tracks from different PVs, two mQv- 
distributions are constructed. One distribution is formed 
from vertices which do not include a so-called large-angle 
track, i.e. the average angle between a constituent track 
and the rest of the tracks in the vertex must be less than 
one radian. The four-momentum of a randomly-selected 
track in the event is added to the vertex to build the second 
distribution. Both of these distributions are taken from a 
large control sample of jet-triggered events. The relative 
normalisations of these distributions are determined using 
an estimate of the "crossing probability" . This is obtained 
by fitting a function representing a Kg mass peak to the 
invariant mass distribution of all two-track combinations 
within three-track vertices, and comparing the integral of 
this function (which is a measure of how many real Ks 
were crossed by a random track to form a three-track ver- 
tex), to the number of two-track Kg mesons observed. 

With the shape of the WDV-distribution determined in 
this manner, the absolute normalisation for the muon- 
triggered sample can be found. It is extracted from a 
maximum likelihood fit, which takes as input the num- 
bers of vertices in each material layer (i.e., the beam pipe 
and three pixel detector layer^ , and the numbers of low- 
track-multiplicity vertices in the regions in which there is 
no material, which are termed "air gaps" . It has been em- 
pirically observed that in each material layer or air gap, 
the relation between the numbers of n-track, (n-f l)-track, 
(n -|- 2)-track vertices etc., is well modelled by an exponen- 
tial function. Furthermore, the ratio of the number of n- 
track vertices in a material layer to the number in the air 
gap immediately outside it, is consistent between several 
different control samples. These relations can therefore 
be used in the fit to predict the number of vertices with 
at least five tracks in the air gaps. This procedure was 
tested on various control samples. For instance, in events 
that fail the muon trigger requirement, the fit predicted 
105. 4± 5.7 events, whereas 110 events were observed. The 



detailed study of the material distribution within the ID is 
provided in Ref. |29| . 
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mDV distribution, scaled by the prediction, is then inte- 
grated to obtain an estimate of the number of vertices with 
TTiDv > 10 GeV, and a final scale factor is applied, based 
on the fraction of events in our muon-triggered sample 
that pass the offlne muon requirements. The scale factors 
are data-derived and arise due to the imperfect simulation 
of the muon triggers; they change the predicted number 
of vertices by up to around 4%. Finally, the background 
from hadronic interactions in the air gaps is estimated to 
be 3.7^3:7 X 10-3 vertices. 

The total background estimate is therefore taken to be 
4^4° X 10"3 vertices. 

7. Systematic corrections and uncertainties 

Several categories of uncertainties and corrections are 
considered. The values of the uncertainties depend on the 
neutralino lifetime and the signal sample under investi- 
gation. However, they can be ordered in approximately 
decreasing size. Uncertainties on the muon reconstruction 
efficiency are largest, followed by tracking and trigger effi- 
ciency uncertainties. The imperfect simulation of pile-up 
and the limited number of events in the signal samples 
lead to uncertainties of a similar size while the integrated 
luminosity makes the smallest contribution to the total 
uncertainty. 

To estimate the uncertainty on the efficiency as a func- 
tion of cr (Fig. [3| from finite MC sample size, the effi- 
ciency in each bin of the two-dimensional efficiency map 
(Fig. |4| was varied randomly within its statistical uncer- 
tainty and the total efficiency was repeatedly re-measured. 
The spread in those results is taken as the contribution 
to the systematic uncertainty on the efficiency at each cr 
value. 

To compensate for the different z-distributions of the 
primary vertices in data and MC simulation, a weight is 
applied to each simulated event such that the reweighted 
zpv distribution matches that in data. This weight is ap- 
plied to both the numerator and denominator in the effi- 
ciency calculation. 

Similarly, a weight is applied to each simulated event 
such that the distribution of (/z), the average number of 
interactions per bunch-crossing, matches that in data. An 
uncertainty associated with this procedure is estimated by 
scaling the {fi) values used as input to this correction cal- 
culation by a factor 0.9 (motivated by taking the differ- 
ence between scaling the MC simulation by (/i) calculated 
from luminosity measurements and by using the number 
of reconstructed PVs in data and MC simulation). The 
difference in the efficiency is evaluated as a function of cr, 
and is applied as a symmetric systematic uncertainty. 

A trigger efficiency correction and its associated system- 
atic uncertainty are derived from a study of Z — > /i"*"/^" 
events in which at least one of the muons was selected with 
a single-muon trigger. The trigger efficiency for selecting 
a signal event is about 90%. The ratio of the trigger effi- 
ciencies in data and simulation is applied as a correction 



factor. Furthermore, the statistical uncertainty on this 
ratio is added in quadrature to the differences in trigger 
efficiency between the Z — >■ jJ,'^ fi^ sample and the signal 
samples to estimate a systematic uncertainty. 

The systematic uncertainty on the track reconstruction 
efficiency is estimated by comparing the track-finding per- 
formance in data and MC simulation. The uncertainty is 
estimated from a study where the number of Ks vertices 
in data and simulation are compared over a range of radii 
and ry-values. Based on the outcome of this study, some 
fraction of tracks in the signal MC are randomly removed 
from the input to the vertexing algorithm, and the same 
procedure is used to obtain efficiency-vs-cr is performed. 
The difference between this and the nominal efficiency at 
each CT point is taken as the systematic uncertainty. 

Since the distribution of the true values of do for cosmic- 
ray muons is flat over the limited do range considered here, 
as predicted by cosmic-ray simulation, the shapes of the 
measured (and simulated) do-spectra of such muons are 
used to determine the accuracy to which the simulation 
reproduces the muon-finding efficiency as a function of do . 
The ratio of the number of observed cosmic-ray muons 
(recorded during 2011) and the simulated muon-finding 
efficiency is studied as a function of do- The degree to 
which the distribution is not flat is used to estimate the 
c?o-dependent muon-finding efficiency, which is used for the 
results in Figs. [3] and |4] The uncertainties on the effi- 
ciency (largely due to the limited statistical precision of 
the cosmic-ray background sample) are propagated to un- 
certainties of between ±3.5% (MH) and ±8% (ML) on the 
signal reconstruction efficiency, depending on the signal 
sample. 

8. Results 

Figure [s] shows the distribution of todv vs. Nj^y for 
vertices in the selected data events, including vertices that 
are rejected by the selection on moy and N^y. The sig- 
nal distribution for the MH sample is also indicated. The 
signal region, corresponding to a minimum number of five 
tracks in a vertex and a minimum vertex mass of 10 GeV, 
is also shown. No events are observed in the signal region. 

Given the lack of candidate events in data, upper limits 
are evaluated on the supersymmetry production cross sec- 
tion (cr) times the square of the branching ratio (BR^) 
for produced squarks to decay via long-lived neutrali- 
nos to muons and quarks. The limits are presented for 
different assumed values of squark and neutralino mass 
and CT, where r is the neutralino lifetime. Based on 
the efficiency, background estimate, and luminosity of 
4.4 fb-^(with the uncertainties on each of these treated 
as nuisance parameters), the CLg method (30) is used 
to set limits on the plane of a ■ BR^ versus cr. The re- 
sulting limits are shown in Fig. ^ with the PROSPINO- 
calculated cross sections for squark masses of 700 GeV and 
1500 GeV. The PROSPINO predictions are shown as a 
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Figure 5: Vertex mass (rrenv) vs. vertex track multiplicity {N^y) 
for DVs in events with no muon requirements other than the trigger, 
and all other selection criteria except the mj^y and N^y require- 
ments. Shaded bins show the distribution for the signal MC MH 
sample (see Table and data are shown as filled ellipses, with the 
area of the ellipse proportional to the number of vertices in the cor- 
responding bin. 



band which represents the variation in predicted cross sec- 
tion when using two different sets of PDFs (CTEQ6.6 and 
MSTW2008NLO [H]) and varying the factorisation and 
renormahsation scales each up and down by a factor of 
two. Since essentially no background is expected and no 
events are observed, the expected and observed limits are 
indistinguishable. Based on the observation of no signal 
events in a data sample of 4.4 fb^^, a 95% confidence-level 
upper limit of 0.68 fb is set on a ■ BR^ multiplied by the 
detector acceptance and the reconstruction efficiency for 
any signal vertex. 

9. Summary and conclusions 

An improved search is presented for new, heavy par- 
ticles that decay at radial distances between 4 mm and 
180 mm from the pp interaction point, in association with 
a high-transverse-momentum muon. This search is based 
on data collected by the ATLAS detector at the LHC at 
a center-of-mass energy of 7 TeV. Fewer than 0.06 back- 
ground events are expected in the data sample of 4.4 fb^^, 
and no events are observed. Limits are derived on the 
product of di-squark production cross section and decay 
chain branching fraction squared, in a SUGRA scenario 
where the lightest neutralino produced in the primary- 
squark decay undergoes an i?-parity violating decay into a 
muon and two quarks. The limits are reported as a func- 



Figure 6: Upper limits at 95% confidence level on a ■ BR^ vs. 
the neutralino lifetime for different combinations of squark and neu- 
tralino masses, based on the observation of zero events satisfying 
all criteria in a 4.4 fb~^ data sample. The shaded areas around 
these curves represent the ztltr uncertainty bands on the expected 
limits. The horizontal lines show the cross sections calculated from 
PROSPINO for squark masses of 700 GeV and 1500 GeV. The shaded 
regions around these lines represent the uncertainties on the cross 
sections obtained from the procedure described in the text. 



tion of the neutralino lifetime and for a range of neutralino 
masses and velocities, which are the factors with greatest 
impact on the reconstruction efficiency. Limits for a vari- 
ety of other models can thus be approximated from these 
results, based on the neutralino mass and velocity distri- 
bution in a given model. 
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